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The Thermal Stresses in Solid and in Hollow Circular Cylinders 

concentrically Heated. 

By Charles H. Lees, D.Sc, F.R.S. 

(Received May 2, 1922.) 

1. The large differences of temperature which exist between the inside and 
outside of the heating or firing tunnels used in pottery and other works, and 
between the inner and outer portions of the concrete or masonry pillars 
supporting the floors of a building in which a fire occurs, make the thermal 
stresses approach the breaking strength of the materials and they are crushed 
or torn asunder. The same may be said to a less degree of chimney stacks, 
and of the walls of the cylinders of internal combustion engines. 

In order to prevent this, it is necessary to design the structure so as to 
withstand the stresses which will be produced by the differences of tempera- 
ture to which they will be subjected. The calculation of these stresses has 
only been carried out for a hollow cylinder of constant elasticity and heat 
conductivity, with the temperature throughout it following the steady 
logarithmic law and the material expanding according to the linear law.* 

It is proposed to treat the problem in a more general way in the present 
communication, and to show that for any temperature distribution, with any 
law of expansion, the stresses may be found by a simple graphical method 
so long as the variation of the elastic constants with temperature may be 
neglected. 

Determination of the Stresses. 

2. Let a hollow circular cylinder, of internal radius r x and external radius r 2 > 

and considerable length, originally at constant temperature, and not subjected 

to any pressure, have its temperature raised so that 

the temperature of each point of a cylindrical 

surface of radius r within it is raised to 0, where 

is a function of r. Let the cylinder then be subjected 

to a uniform pressure, Ei, at its inner, a uniform 

pressure, Rs, at its outer curved surface, and a 

uniform pressure, X, at each of its end surfaces. To 

find the mechanical and thermal stresses at each FlG * L - Radii > r i and r 2 > 

sectional areas, s t and s 2 . 

point of the material of the cylinder. 

In the first instance we note that a uniform rise of temperature would 
produce a uniform expansion of the cylinder, and no stresses would result 

* See R. Lorenz, * Zeitsch. Yerein Deutsch. Ing.,' vol. 51, p. 743 (1907). 
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from such an expansion. We therefore divide the final temperature, 0, into 
two parts : first, a uniform temperature, 0, to which the cylinder would have 
to be heated to give the same expansion of the whole shell as that produced 
T)y the actual temperature, ; second, the excess, 0', of the actual tempera- 
ture over this uniform temperature. 

Thus = + 0'. 

Whatever the character of the law of expansion of the material of the 
shell with rise of temperature, we may write v = Vo(l+ot e ) where v is the 
"volume of an element at 0° C. and v e at 0° C, and a Q is the increase of volume 
of 1 c.c. when raised from 0° C. to 0° C. Hence if a is the mean expansion for 
the whole shell and a/ the excess at a point where the temperature is over 
the mean expansion, we have ot = ar-f a/. 

Next we notice that, although the pressures on the end surfaces are 
"uniform the pressures in planes parallel to them cease to be uniform, owing 
to the variation of 0' with r. In planes at distances from the end surfaces 
greater than the thickness of the cylinder these pressures become indepen- 
dent of the positions of the planes. 

In the middle length of the cylinder points initially in a plane perpen- 
dicular to the axis of the cylinder at x are still in a plane perpendicular 
to the axis and at x+x' after the expansion, and the initial distance, dx 
apart of two such planes becomes (l+(dx/dx)dx after expansion, and is 
independent of r. The same statement would hold for the whole length of 
the cylinder, provided the end surfaces were kept plane by pressures or 
tractions applied in a suitable way over them. 

Under these circumstances the only stresses acting on an elementary 
portion of the material bounded by two planes perpendicular to the axis at 
distance x and x + dx from some fixed plane, by two cylindrical surfaces of radii 
r and r+dr, and two planes through the axis making angles (j> and <£-f c£<£ 
with some fixed' axial plane, are the equal axial pressures, X, acting on the 
surfaces perpendicular to the axis, the radial pressures, R and R-f dR, on the 



c 



R 




It+dft 



Fig. 2.— The pressures- 



-R, radial ; T, tangential ; and X, axial, at the faces of an 
element. 
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inner and outer cylindrical surfaces respectively, and the equal tangential 
pressures, T, acting on the axial surfaces at right angles to those surfaces 
(fig. 2). If the displacement outwards of a point originally r from the axis 
is r\ the stress-strain equations for the elementary portion of the material are 

I |_T + <r(E + X)| =^-W, (2.2) 

K _x+o ' (E+T) } = £" ia/ ' (2 - 3) 

where e is Young's modulus, a Poisson's ratio for the material, and dx r jdx is 
independent of r and x and may be written a x . 
The stress equilibrium equation is 

T = K + r^P, (2.4) 

Multiplying (2.3) through by a and adding to (2.2) we get 

I j"-T + o-K-f<r 2 (E + T)T =£ +<ra x -l£? . a,'. (2.5) 

e ' L J r o 

Subtracting (2.2) from (2.1) we get 

l + <r (T _ R) dr' r' 



e * dr r 

Substituting for T—E from (2.4) this gives 

1 + cr dB, _ d it' 

e dr dr\r_ 

and ~ R = _ + a constant. 

e r 



(2.6) 



Differentiating (2.5) with respect to r, assuming that the variations of 
€ and <r with the radius are negligible and subtracting the result from (2.6) 
we get 

1 — <r 2 d ,m , ipx ___ l + <rdoLe' 
€ dr 3 dr ' 

Integrating this gives us 

T + E = | r ^-a/ + A / , (2.7) 

6 1 — G" 

where A' is a constant of integration. 
Substituting for T from (2.4) we get 

2R + r^ = £-*-*,' + A'. 
dr 3 1— ar 
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Multiplying through by r and integrating we have 

t*R = B' + A'~ + 1 r-2- ( a / . rdr. 

Jj OX — <T J 

The form of this equation suggests the expression of the stresses in terms 
of § the area of the circle of radius r instead of r. The equation becomes 

Es = As + B + g j^- Pa/ . ds, (2.8) 

ttA' 
where A = — — and B = ttB'. 

Hence Ei and E 2 , the pressures on the inside and outside surfaces 
respectively, satisfy the equations 

1 e f 5] 
Ei$i = Asi + B + g- r~— oce ds, 

1 e f s 2 
E 2 s 2 = As 2 + B + ^ . oce ds. 

O 1 — or 



1 



olq ds 



l e 
Hence Es—EiSi = A(s— *i)+« 7 . 

O x — - O" 

and E 2 s 2 --RiSi = A(s 2 — s x ) since a/ds = 0. 

On elimination of A these give 

p _ j^ggg — E1S1 (Ei — E 2 ) sis 2 1 _, 1 € E' , . 

s 2 — $i 52— «i s 6 1 — a s 3 

where 

E'= L;js= (a* — a)ds = E — £($— Si) = E — E 2 ±ZfL 

J«, J«i ' S 2 — Si 

E 2 being the value of E at the outer surface. 

Or if u be written for 1/s throughout the coefficients 

P — ^2^1 -— Ei U2 , Ei — E 2 1,1 e E' ,q 1 n\ 

wi-— U2 Wi— ^ 2 s 6 1 — 0- s 

^E 
Equation (2.4) becomes T = E + 2s_ hence 

T = ^^i-Ei^^Bj-Ej 11 eU,_W\ (2J1) 

U\-~ U 2 U1 — U2 s 1— cr\ sy 

On substituting for E and T in equation (2.3) we get 

X = - € a x + 2<r' R2Ul -' RlU2 + l r ^ 2a*'. (2.12) 

U\—U2 O 1 — cr 
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On substituting for E, T, and X in equation (2.2) we get 

r 2 s € L u\—v* tii—u 2 sj 61-o-s 

(2.13) 

The fractional increase of length a x in the axial direction and the part X x 
of X which is independent of r and x may be given any arbitrary values so 

long as their sum X x + ea x is kept equal to 2cr 2 x ~~ l 2 . 

U1 — U2 

Thus if the increase of length a x is taken zero X x = 2a — — — , while if 

Ui—'ii 2 

X x is taken zero a x = 2a - 2 1 "" — i-2 

Whatever value is assigned to a x it only affects the expressions for the 
mechanical portions of the stresses and displacements, the terms for the 
thermal stresses and displacements are unaffected. They depend on E' and 
dW/cls = a e ' only. 

The mean pressure 4-(B + T + X) at the point r is given by 

I (E + T + X) = % (1+cr) Jh&=Ml~l a x +\<**L t (2.14) 

the shear stress in the plane T, X is given by 

1 /y t\ 1 — 2cr E 2 ^i— Ei^ 2 , 1 Ei— B 2 1 e rt ,1 e IE' 

3 b Ui — u 2 3 ux—*u 2 s 3 6 1— a 3 s 

(2.15) 
and in the plane E, T by 

1 /T) T\ ^ 1^1 — -^2 1.1 € / 2 / 2 E \ /oi/»\ 

3 3 ^1—^2 5 6 1 — <r \ 3 3s/ v . 

and in the plane E, X by 

■ 1 /"D Y\ 1 — 2<T E.2% — Ei1^2 , IE1-E2 1 , 6 

o \**— <&-) n — • r^ ■ . - + 77 $z 

o 6 U1 — U2 3 Ux—u 2 s 3 



+ -Jrb(iT-i-'} (2 - 17 > 

The terms involving the internal and external pressures are identical with 
those given by Lame in his Leqons and quoted by Love. 

If the cylinder is solid, the constant of integration B in equation (2.8) 
becomes zero. On making this change in the calculations giving the stresses 

it is found that the term 2-LZ lU2 in the expression for each stress reduces 

u x — u 2 

to B 2 — the pressure on the outside surface — the term involving Ei~B a 

* Love, * Theory of Elasticity,' 2nd edit., p. 141 (1906). 
VOL. CI.— A. 2 E 
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disappears, that involving a x remains, while the terms giving the thermal part 

of the stresses remain unchanged in form, but E' is now I a/ * ds and — at 

Jo 5 

that is <%/ at the centre, so that the thermal radial 



s = becomes 



ds /o 



pressure does not become zero at the centre, but the shear stress in the E, T 
plane does. 

Graphical Constructions for the Stresses and Displacements. 

3. The foregoing expressions for the stresses and displacements lead to the 
following methods of representing and determining both graphically. 

A. The Stresses Due to the Pressures on the Inner and Outer Surfaces. 

Set off along the axis of abscissae, fig. 3, the reciprocals OTJi = u x and OU2 = u% 
of the areas s\ and s 2 of the circles bounding the normal section of the 
cylinder. At IT 2 erect the ordinate U 2 E 2) equal on some convenient scale to 
the pressure E2 on the outside surface. Join the end E 2 of this ordinate to* 
the point Ui, and produce to cut the axis of pressure in So. Through So 
draw a straight line parallel to the axis of abscissae, and let it cut the ordinates 
through U 2 and Ui in S 2 and Si. Set off above Si a length S1E1, equal to Ei 
the pressure on the inner surface, join the end of this ordinate to S 2 , and 
produce to cut the axis of pressure in So'. 

Then OS is E 2 . — ^— , SqSo' is E x . J02-, hence OSo' = ?^LZ^1^, that 

Ui — U2 U\ — U2 U\ ~~ U2 

is the constant term in the expressions for the stresses. 

Ei — E 2 



To get the term involving 



u 7 erect at Ui, fig. 4, the ordinate 



ui—u 2 
UiA = Ei — E 2 , join the end of it to the point U 2 and draw through a 

straight line OB parallel to this join. The ordinate US at any point U is. 
then = 5LZ^£ 



u. 



ui — ^2 





—d B-r^ 



Fig. 3. 
Fig. 3. — 0U 2 = %, OUj = u l9 OS ' the mean pressure (tension) determined by 

Fig. 4. — The variable stress US determined from UjA = R x -R 2 , U 2 A and OSB. 
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The mechanical stresses are then given as follows : 

E = 080' + US, 

T = OSo'-US, 

X = — ea x + 2<r . OSo', 

T € 

the mean pressure 



{-(l-2o-)OS , + US}-(7^, 



the shear stresses 



and 



I (K + T + X) = |(l+<r) OS '-| a x , 



\ (X-T) = -^~ OSo'+i US- I « 2) 

i (K-T) = I US, 
I (R-X) = -L^ OSo'-^ US + I a x . 



(8.1) 

(3.2) 

(3.3) 
(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



B. Stresses Due to the Hxcess of the Dilatation at a Point over the Mean 

Dilatation. 

Set off as abscissse OAi OA2 the areas s\ and s 2 of the inner and outer 
boundaries of the section of the cylinder. As ordinates take the temperature 9 
of each narrow ring of internal area s, and draw the curve of temperature 
distribution throughout the section, fig. 5. 

For each temperature estimate from tables of dilatation of the material, the 
dilatation &$ of 1 c.c. of each part of the section, and with a e as ordinates, draw 
the curve of dilatation per cubic centimetre D to correspond to the 
temperature curve. 




Fig. 5. 





Fig. 7. 



Fig. 5. — OAj = s ly OA 2 = s 2 , 6 the curve of temperature, and a Q that of dilatation. 
Fig. 6. — The expansion curve E, OA x = s 19 OA 2 = s 2 . 

Fig. 7..— AjE'Aa the excess expansion curve ; OC parallel to the tangent at P ; 

ADC a convenient ordinate. 

2 E 2 
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Find the total expansion E of the shells within the radius r, if they had 
expanded independently of each other, that is, find the area enclosed between 
the axis of area, the curve- D, the ordinate through Ai and that through s for 
different values of s up to s 2 , and plot the values obtained as ordinates against 
the areas as abscissae. Join the end of the curve at E2 with its starting point 
at Ai, fig. 6. 

The excess of the ordinate of the curve at any point over the ordinate of 
the join is the ordinate of the excess expansion curve W required, fig. 7, the 
scale of ordinates being for convenience increased. 

Join P a point on the curve to the origin 0. Through draw OB parallel 
to the tangent to the E' curve at P and at A, corresponding to area s, a 
convenient multiple of 10, erect a perpendicular to the area axis cutting OB 
in C and OP in Di. 

Then the thermal stresses are given as follows : 

E - J « . AD/OA, (3.1') 

T ■= 1--1-. (2AC-AD)/OA, (3.20 

O 1 — cr 

X = I --L. 2 A£, (3.3') 



6 1-a OA' 
r 6 1-<rOA' 



(3.40 



the; mean pressure 



l (R+T + X) = ^4.^ (3.50 

the shear stresses 

1 , v „. 1 ■ e 1 AD 1 /p T , _ 1 e 2 AD -AC 

(3.6' & 1') 

It will be noted that for the thermal stresses E -f T = X at any point, and 
since at the surfaces E = 0, T = X at them both. 



Numerical Evaluation of Stresses in Practical Cases. 

4. In applying the foregoing theory to practical problems it seems desirable 
to select in the first instance cases of a comparatively simple type in which 
the stresses due to mechanical causes are small. Two such cases are (1) that 
of a long pottery kiln, a section of which is kept at the firing temperature 
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either permanently or alternately with its neighbours, and (2) that of a 
circular pillar of concrete or of masonry supporting the ceiling of a warehouse 
in which a fire occurs. 

Case (1) along Pottery Kiln or Furnace. 

5. The distribution of temperature throughout a cylinder concentrically 
heated and of length sufficient to make the effect of the ends negligible is 
determined by the equation 

w = lM rh %A\ (51) 

dt r dr\ drj 

If, in the first instance, we take the density p f the specific heat c and the 
heat conductivity k as constant, and the temperature distribution throughout 
the wall of the furnace as steady, 0\ being the temperature of the inner 
surface at r h 2 that of the outer surface at r 2) we have 

0—02 = logjv-logg = logs 2 — logs ,gg\ 

0i — 02 logr 2 — logn logsa — logsi' 

For simplicity we may take the dilatation a e proportional to the rise of 
temperature, that is a = a0. The expansion E of the rings between 
§i and s, if each expands independently of its neighbours, is given by 

E = a(0~0 2 ) ds, and may be found either by graphic integration of the 
curve of temperature or directly from the equation 

E = a (0i— 2 ). r— {s (log sa — logs) — si(logS2--logSi)-M— §i}, (5.3) 

logs 2 — logsi 

and plotted as in fig. 8, where the temperature and the coefficient of a(#i— #a) 
in E are shown for a cylinder when si = # 75s2. In that case the temperature 
curve is nearly linear. 

The total expansion of the area of section of the kiln wall is given by the 
ordinate at A 2 of the E curve, fig. 8, and has the value 



logs 2 Si L Si J 



s — Si 



The excess expansion E' = E— E 2 is represented by the excess of the 

S2 — Si 

ordinate of the E curve over that of the join of the two ends of that curve in 
fig. 8 and is given by 

t?' /z) a \ f logs 2 — logs n s 2 ~ si /K , N 

E = a (01—02) < s * ' & — si.. y . (5.4) 
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The excess of the dilatation at any point over the mean dilatation for the 
whole wall is represented by the slope of the E' curve, and is given by 

as L logss— logsi S2—S1 

On drawing the lines OP, OB and the vertical ADC, corresponding to fig. 7 5 



10 for E 




-1-0 
Fig. 8. — Temperature^ expansion and stress curves for a hollow cylinder heated within. 

or on substituting the values found above in the expressions for the thermal 
stresses and displacements, we get 

1 e AD _ 1 e / gi _ ft) f logs 2 -~logs __ s 2 /s-l \ ^ 55) 

6 1 — cr Llog$2— logSi s 2 fsi--lj 



E 



6 1 



OA 



Thus the radial pressure is zero at the inner surface, rises to a maximum 



6 1 



.a(^-ft) 



-t-^—1 — 



(S2/Si)-~l S X log S 2 — log Si 



= I _1_ « ld 1 -0 2 ) f S2 / s + 1 2-(Iog* a -logg) "\ ^ 
6I — 0- LS2/ 5 i — 1 log $2— log Si J 9 
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at the section s = sis% S 2 ~"" — s_i ? or s = s* , f ?' V > an d becomes zero at the 

S2—S1 (S2/S1)— 1 

outer surface. If the outer temperature # 2 is higher than the inner &i these 

pressures become tensions. 

In the same way the circumferential pressure T is given by 

6 1 — o- OA 6 1 — <7~" v " x "*' I.S2/-1 - — & 

(5.6) 
Thus the circumferential pressure is 

1 6 „(& ft \ I 2s2 2 

1 — a l$2 — Si logs 2 — log^i 

at the inner surface T\ t becomes zero slightly within the ring of sectional area 

($i-f-$2)/2,andis— -—-i~ a(#i — 6 2 )\ -r- — —+- >, that is a tension 

' 6 1— <r L S2—S1 logs 2 — logSiJ 

at the outer surface r 2 * 

The longitudinal pressure X is given by 

6 1 — c OA 6 1 — or Llog s 2 /si s 2 —SiJ 



It changes from - -~^ — a (0i — #2) 2 -< — r > at the inner surface 

6 61-0- V1 J/ Aa-*i logs 2 /sJ 

to zero at the surface — (-2) and becomes 

6 1 — <r Ll0gS2/5i 5 2 — Si J 

at the outer surface. 

The thermal displacement other than that due to the mean rise of 
temperature is given by 

!l 1£ 1^^ l^. tt ^ (5.8) 

r 2 s 6 1-<rOA 6 1-<r v Uog« a /«i s 2 /si — lj V ; 

Thus the excess thermal displacement is zero at the two surfaces and has 

its maximum value -- — - . ol(&\— 2 ) , r~ \ log ( , * . ) — 1 > at the 

6 1 - a v 'log s 2 /si I & \log s 2 Si 1 J 



surface - = 7 * z . 

S2 S2/S1 — 1 



To this displacement must be added that due to the mean rise of 

temperature which is given by — - = - — = - a (#i — $2) < g2 "" gl ^s x 
* b J r 2 $ 3 - J Llog 82/81 

* The values of the tangential and axial stresses at the two surfaces are identical 
with those given by E. Lorenz, loc. cit. 
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The coefficients of - — — . a {61 — 62) in the expressions for the stresses 

are shown in fig. 8 as functions of s for the case s x = 0*75s2- If 

61 — 62 = 1000° C, the coefficient of cubical expansion of the material 

a = 12x10" 6 , Young's modulus e = 10 11 and Poisson's ratio or = -, the 

o 

value of -- ol {61 — 62) is 3x 10 8 dynes per sq. cm. The tensile strength 

o 1 — (J 

of fire brick is of the order 2xl0 7 dynes' per square centimetre (300 lbs. 
per square inch). The tensions T and X in the outer portion of the wall 
are, therefore, capable of producing cracks in the material which will be 
radial, and, since at the outer surface T and X are equal, will extend in the 
longitudinal as often as in the transverse direction. 

6. If the temperature -throughout the section is not steady but variable, the 
number of possible distributions is infinite, but of these cases two are worth 
investigating as of practical value (A) when the kiln or furnace initially at 

62 has its inner surface raised at a given instant to 6%, and (B) when a kiln or 
furnace has been in operation with its inner surface at #1 and its outer at 62 
for a considerable time, and at a given time has the temperature of its 
inner surface reduced to 62. 

A. The kiln is brought into service. 

For the first few minutes the only change of temperature is that close to 
the inner surface. The mean temperature is scarcely changed from 62, the 
E' curve has very small ordinates but 6' has approximately the value 61 — 62 
close to the inner surface. 

Hence at this stage, close to the inner surface, E the radial pressure and 
v' \r the radial displacement are very small, but T the tangential and X the 
longitudinal pressures are both very nearly double what they finally become 
in the steady state, that is near the inner surface they both have initially 
the value 

and decrease to half that value as the temperature distribution approaches 
that of the steady state. 

B. The kiln is put out of service. 

For the first few minutes the only change of temperature will be at the 
inner surface. The mean temperature of the section is scarcely changed and 
the E' curve has only a very smail negative ordinate near its zero value at 
the inner surface. The approximate value of 6' at the inner surface is 
changed from (6 1 -6 2 )/2 to -{6 l -6 2 )/2. 
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Thus E the radial pressure and r' jr the radial displacement are very small 
close to the inner surface ; but T the tangential pressure and X the ; 

longitudinal pressure change from the approximate value - ■ ■ (61-^62) to- 
ol — or 

I G.OL 

{61 — 62), that is from pressures to equal tensions. 



6 l-<7 

We have seen that, for a furnace with a difference of temperature 
of 1000° C. between inside and outside surfaces of the walls, composed? 
of material with a coefficient of cubical expansion a . == ■ 12 * l-0~* 6 j.' a 
Young's modulus € of 10 11 and a Poisson's ratio a of ^-, the value of 



€U 



(61 — 62) is 1'5 x 10 8 dynes per square centimetre (2250 lb. per square 

12 1 — cr 

inch), w 7 hile the tensile strength of the material of the walls is of the order of 
2x 10 7 dynes per square centimetre (300 lb. per square inch),, and that of the 
mortar about one-fifth of this. The sudden cooling of the interior of a 
furnace such as that described above will therefore subject the interior of 
the wall to tensile stresses in both longitudinal and tangential directions 
much above the tensile strength of the material, and radial cracks will be 
produced at the inner surface, which, since the stresses are equal in the two 
directions, will run indiscriminately longitudinally and tangentially. 

Through the kindness of Sir Eobert Hadfield, I have been able t6 inspect, 
the walls of furnaces which have been in use for various periods at the 
Hecla Works> Sheffield, and I find that the above results are confirmed by 
experience. 

Case 2. — A Concrete or Masonry Pillar subjected to Fire. 

7. A circular pillar of radius a, initially at a uniform temperature 0°, has> 
its outer surface suddenly raised to 6 a . 

Assuming that the physical constants of the material of the pillar may be- 
taken as uniform throughout the material, the distribution of tempera- 
ture throughout the centre length of the pillar at any subsequent time will 
be given by the equation 

M-lll(r^L) (71> 

dt cp r dr \ dr/ 

The solution of which under the conditions laid down is 

6 = e a - 26 a 2 e~*i"W"r* . Jo Jfyi2\ * (7.2> 

where Jo (|S) is the Bessel function of zero order, Ji (0) that of the first 
order, and the /5's are the roots of the equation Jo(/3) = 0. 

* Compare Byerley, ' Fourier Series and Spherical Harmonics,' p. 2291 
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For a disc of the pillar initially 1 cm. thick cut perpendicular to the axis 
we therefore have 

JO L^ ™n JlKPn) J 



Hence at r = # 

A 



and 

E' = E-E a . ^ = ^eate-MM** . Jt fi^iiMa) ^ 5 



n 



a 3 A, 2 I /3»r/a Ji(/?„) 

Hence the stresses are given by the equations 
E = \ -_L_ . E. = | . _i_ . *e a 2r*K*U»t 

6 l — O- 7TT J 3 1 — 0- 

_L J" 1 _ 'h(Pn r M &» 

T = i — i-./2 -^ — ^ = - 6 a^Se-WW *'' 

6 1 — <7\ 27rrrfr 7rr 2 / 3 1 — <r 



1 



A 2 Ji08») 



./3 n \ a J r \ a 



1 € c\ dht Z € 



o 1 — cr iirrar 6 1 — a 



/3 n 2 /3 n Jl (Bn) W 

and the expansion by the equation 

r ' — * s _ 2 1 + cr ^ aQ^-klcpWa) t 1 



(7.3 



\ 



(7.4) 



(7.5) 



Jo(^r) k (7.6) 



r 2 5 3 1-cr A 2 

^lifinrja) ^ .J& \ + 1 a a |l_42 i ^»1, (7.7) 
/3 n f/a Ji(/3 n )J 3 L Pn J 

the last term representing the expansion due to the mean rise of temperature. 

It will be noticed that T-fK = X, that is, the stress in the axial direction at 
a point is equal to the sum of the radial and tangential stresses at the point. 

8. These expressions are convenient solutions of the problem when t, the 
time which has elapsed since the outer surface of the pillar was first heated, 
is sufficiently large to reduce each series to its first term, but, in general, the 
series converge very slowly, and many terms must be taken before an 
accurate numerical result is obtained. 

As it is in the earlier stages of the heating, before the temperature near 
the middle of the pillar has commenced to rise, that the thermal stresses 
appear to produce their greatest effect, it is desirable to attempt to represent 



Hollow Circular Cylinders concentrically Heated. 425 

the temperature distribution for small values of t by a simpler expression, 
which will reproduce the principal characteristics of the actual case. 
We have found 

@n 2 L #»r/a Ji (>8») J 



*W / 2a 1 2\2a/^12\2a/ 



Substituting we find 
E' = 47ra(9 a S6-^(^) 2i 

** -iQgn/2) 2 (a 2 ^^)/^ + ^(/3 n /2)H^-^)M 4 -». (81) 
*/3* 2 l-i(/3«/2) 2 + T V(ySn/2) 4 ' 

and on substituting s/s a for T^/a 2 

W — — nfi c y P -klcp • (j8,/a)2i ^ ( s a~' s )/ S a~-Tt (@n / 2)' (Sg -**^) / Sg — ... 

°- 2 • i-h/w+tVOW—. • 

The form of this expression suggests the assumption that at any time 
during the early stages of the heating process the value of E may be taken 



oc0 a I s\ n 



as approximating to — ^ . 5 (— i where n has a value which decreases as t 
increases. Then 



E' = -- 2fi-«. 

% + l 



A/ J 



(8.2) 



0' = -^~(l-(n + l)(±)"\, d= ° a 



%+l [ \sj J n-hV 

and the values of the stresses and displacements are readily seen to be 

B * « ^jl-fifj, (8.3) 

6 1-<r w + 1 I W J 

T Hrira{ 1 -< 2 » + <)*}' < s - 4 > 

the last term giving r'/r for the mean rise of temperature &, 

Thus the mean temperature B increases as n decreases. 

The radial stress E is zero at the surface, is a tension within, which 

1 e ad 

increases in value rapidly towards its maximum value—- . — ~ at the 

b 1 — <r %-h 1 



426 Dr. C. fit. Lees. The Thermal Stresses in Solid and. in 

centre and has nearly that value for a considerable distance from the centre 
if n is large. 

The tangential stress T has the value - . — ~ 2n and is a pressure 

& 6 l-o- , n + 1 r 

I s \ n 1 
at the surface, decreases to zero at the section ( - = ~ becomes a tension 

\sj 2n + l 

1 e ol6 

and has the value—--— — . . — ~ at the centre, attaining nearly that value 

6 1 — a n + 1 ' 

at a considerable distance from the centre if n is large. 

The axial stress X has the value - e — £- . 2n and is a pressure at the 

o 1— a n + 1 

surface, decreases to- zero at the section ( ~"\ = -- becomes a tension and 

\sj n+1 



has the value — -- — — . — —. . 2 at the centre, attaining nearly that value 

6 1 — .<r n+1 ■ J 

at a considerable distance from the centre if n is large. 

The displacement term r' lr has the value -2-4 /-+•? ' ) at the surface 

n + 1 \3 6 1 — crl 

and the displacement r f decreases to zero at the centre. 

The curves of fig. 9 give the coefficients of - . - — ~ in the expressions- 

: 6 l — o- n + 1 

for the stresses for the case n = 8. The curve for the axial stress X shows 

also the distribution of temperature above and below the mean temperature, 

and with the coefficient ~ - — . — ~ that part of the expression for the 

3 1 — 0- n + 1 k r 

displacement term r' jr which depends on the variation of the temperature 

from the mean temperature. 

During the early stages of the heating, the longitudinal tension X in the 
central portion of the section is double the radial tension E and the 
tangential tension T. If the pillar is not subject to other mechanical 
stresses, and the magnitude of these thermal stresses approaches the tensile 
strength of the material, cracks in planes perpendicular to the axis of the 
cylinder will be more frequent in the central portions than cracks in other 
directions. 

As, however, the tensile strength of most materials is less at high 
than at low temperatures, it is possible that the radial tension, which 
persists, but with diminished magnitude in the outer and hotter layers near 
the surface, where the tangential and longitudinal stresses have become 
pressures, may still be large enough in these layers to exceed the tensile 
strength of the material. In any layers in which this is the case, cracks wiE 
be formed parallel to the outer surface of the cylinder. 
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The tangential tension T, like the radial tension B, is half the longitudinal 
tension X in the central portions of the section, but becomes zero, and then a 



3=0 




Fig. 9. — Temperature, expansion and stress curves for circular pillar, heated outside, 

in case n — 8. 

pressure in the outer portions, where K is still a tension. Hence the tendency 
for cracks parallel to the axis to form during the heating process is less than 
for cracks perpendicular to the axis in the central portions, and less than for 
cracks parallel to the outer surface in the outer portions. 

A Heated Solid Cylinder has its Surface Suddenly Cooled. 

9. In this case the stresses and displacements are identical in magnitudes 
with those determined for the heating of the cylinder, but the sign of each is 
changed, a pressure becoming a tension and a tension a pressure. In fig. 9 
we may consider the signs of the ordinates reversed. At the stage of the 
cooling shown, which we may regard as representative for the earlier stages, 
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the greatest tensile stress is the tangential one, T, which, except just at the 
surface, is slightly larger than the axial tensile stress X, and decreases 
rapidly in magnitude as we approach the centre. While T and X are both 
tensions, the radial pressure E is positive — that is a pressure. During 
cooling, we see, then, that there will be a slightly greater tendency for T ta 
form cracks running parallel to the axis than for X to form cracks running 
round the cylinder. 

There will be no tendency for E to form cracks running parallel to the 
outer surface of the cylinder. 

In comparing the stresses produced on cooling with those due to heating, 
it must not be assumed in general that the contraction on cooling is 
identical with the dilatation on heating. As a rule, building materials 
expand more on first being heated than they do at any subsequent heating,., 
and after each heating contract on cooling by a less amount than they 
expanded. Each heating and cooling results in a permanent increase of 
volume, which decreases in magnitude with each successive change.* In 
such cases it is therefore necessary to use the general theory of Section 2. 

10. It is difficult to find experimental observations with which to compare 
the above deductions. The only tests of cylinders subjected to fire without 
mechanical loading appear to be those of W. A. Hull of the Bureau of 
Standards,! which were carried out on cylinders of clay, concretes, gypsums 
and mortars 8 inches in diameter and 16 inches long, heated in a gas furnace 
for 3 J hours to 927° C. (1700° F.). These cylinders are too short to justify 
the neglect of the effect of the end surfaces, but the temperatures are of the 
order of those used in the numerical work of this paper. Mr. Hull draws no 
general conclusions as to the tendency of the cracks, but the illustrations 
show a tendency for the major cracks to run in planes perpendicular to 

the axis. 

For longer vertical cylinders under mechanical stress due to a load carried 
by the cylinder, the observations of HullJ appear to be the only ones 
available for comparison. These columns were 12 or 18 inches in diameter 
8 feet 9 inches long, of various materials, and were heated for 4 hours in a gas 
furnace to 1100° C. while under load. Both plain and reinforced concrete 
columns were tested. 

In this case in addition to the longitudinal thermal stress X there was the 

* See Day, Sosman and Hostetter, ' Amer. Jour. Sci.,' vol. 37, p. 1 (1914) ; and Mellor 
and Campbell, ' Trans. Ceramic Soc.,' vol. 15, p. 80 (1915), and references given by them. 

t W. A. Hull, "Heat Insulating Properties of Materials," 'Proceedings of the 
American Society for Testing Materials,' Philadelphia, June, 1917. 

% See W. A. Hull, " Fire Tests of Concrete Columns," ' Proceedings of the American* 
Concrete Institute,'* vol. 15, p. 1 (1918). 
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stress due to the load which, in the case of the plain concrete columns, was of 
the order of 800 lbs. per square inch. 

According to the curves of temperature distribution given by Hull, the 
curves of fig. 9 correspond approximately to the state half an hour after the 
heating commenced. With the constants for concrete taken as -a = 4x 10"" 5 , 

1 € OC0 

e = 1 x 10 n , a = 1/3, the value of the coefficient - . — ~ of the equa- 

o 1— a n+1 

tions (8.4, etc.) and of the ordinates of the curves of fig. 9 is 1*2 x 10 8 dynes 

per square centimetre (1800 lbs. per square inch). Thus in the central 

portion of the section of the column the thermal stresses K and T are tensions 

of 1*2 xlO 8 dynes per square centimetre and X a tension of double that 

amount. The latter has to be diminished by the pressure due to the load and 

we get a resultant tension of about 2800 lbs. per square inch, which is greater 

than the radial and tangential tensions. There will, therefore, be a greater 

tendency for cracks to form perpendicular to the axis in the central parts of 

a column than in other directions. In the outer portion of the section the 

tangential stress T becomes a pressure before the radial E, and neither are 

affected by the loading. In these portions E is more likely to form cracks 

parallel to the outer surface than T is to form cracks parallel to the axis* 

The positions of the cracks due to E will be determined by the decrease of 

the tensile strength of the concrete as the temperature rises. The curves of 

fig. 9 show that from s/s a = # 6 to 0*8 while E decreases little the temperature 

increases considerably and the tensile strength therefore decreases. 

Hull's observations and photographs all indicate that the materials crack 
along a surface about 1J to 2 inches within the outer surface (s/s a , 0*6), and 
that the outer shell may, or may not, fall away. The existence of these 
cracks within columns which have been subjected to fire, but which appear 
on the surface uninjured, has been known for some time.* 

As the tensile strength of concrete, at ordinary temperatures, is of the order 
1 x 10 8 dynes per square centimetre (1,500 lbs. per square inch), it is seen 
that the thermal stresses in the heated columns, at the stage of the heating 
represented by fig. 9, are sufficient to account for the cracks produced. 

As the heating continues the tensions in the central portion of the section 
at first increase, while the area subjected to them decreases. When the 
temperature at the centre begins to rise the stresses decrease. 

Where the column is reinforced with vertical rods, with or without 

* See F. J. T. Stewart, * Report on Edison Works Fire,' Scientific and Industrial 
Research Department, British Fire Prevention Committee's Red Book No. 204, p. 28- 
(1916) ; and A. J. and T. D. Mylrea, * Report on Quaker Oats Works Fire,' British Fire 
Prevention Committee's Red Book No. 225, p. 56 (1918). 
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horizontal binders, the crack produced follows the rods or binders, indicating 
that the tensile strength of attachment of metal to concrete is less than that 
of the concrete itself at the temperature attained. 

So soon as the outer shell is detached from the inner portion it has to 
<carry its share of the load without lateral support from the inside. It 
therefore buckles like an overloaded hollow column and eventually spalls off. 

During the cooling process which, in Hull's experiments, was a gradual one, 
;and began long before the column had wholly acquired the higher temperature, 
vertical cracks, starting at the outside surface and extending inwards, would 
be* most likely to form. The tendency would be the greater the more 
rapid the cooling, but less than it would be if the whole column had acquired 
the higher temperature before cooling. The photographs of the columns, after 
the tests given by Hull, nearly all show these vertical cracks,* but there is no 
indication as to whether they were produced during the cooling or not. 

So far as the experimental evidence available goes it appears to support the 
theoretical deductions of this paper, but it is desirable that the thermal 
stresses should be computed either graphically or otherwise in all cases in 
which the necessary observations are available, in order to determine to what 
: extent the explanation of the behaviour of furnace walls and columns here 
offered is adequate. 

* W. A. Hull, * Tests of Concrete Columns, 5 p, 15. An important contribution to our 
knowledge of the effects of fire and water on loaded columns has just been made by 
JMessrs. Ingberg, Griffin, Kobinson and Wilson. See * Bureau of Standards Technological 
Paper No. 184 ? (1921). 



